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Following the great success of the current Imaging Atmospheric Cherenkov Telescopes, the next 
generation gamma-ray telescope arrays, Cherenkov Telescope Array (CTA), is being prepared. 
CTA will have an order of magnitude higher sensitivity and wider energy range than current 
instruments. Large samples of very high energy (VHE; > 30 GeV) sources will be obtained with 
CTA. Here we discuss potential AGN population studies for the CTA era concerning the overall 
source counts and the statistics of high redshift sources. Based on our latest blazar gamma-ray 
luminosity function (GLF) model, we find that CTA will detect ~ 50 and ^160 blazars with 1 
year and 10 years of blank field sky survey, respectively. CTA is also expected to find a blazar at 
z ~ 1.4 (~ 20 blazars above z = 1) based on the blazar GLF. Furthermore, we also examine the 
detectability of high redshift Fermi blazars. By extrapolating the Fermi blazars' spectra with a 
power-law, we find that CTA has a capability to detect a Fermi blazar at z = 2.49. 
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1. Introduction 

Following the great success of the current Imaging Atmospheric Cherenkov Telescopes (IACTs), 
the next generation gamma-ray telescope arrays, Cherenkov Telescope Array (CTA), is being pre- 
pared It is designed to have an order of magnitude higher sensitivity and wider energy range 
compared to arrays. 

Current IACTs have already found ~100 very high energy (VHE) sources, including more 
than 40 blazars up to redshift z = 0.536 1 . Blazars, a class of active galactic nuclei (AGNs), are the 
dominant population in the extragalactic gamma-ray sky. Almost all of the extragalactic sources 
detected by EGRET and Fermi are blazars [§, ||]. In the CTA era, it is naturally expected that 
the number of blazars detected in the VHE range increases dramatically and that a much deeper 
VHE universe will be explored. Therefore, it would be possible to carry out a population study 
of VHE blazars with CTA, which would provide a crucial key to understanding AGN populations, 
high-energy phenomena around supermassive black holes in AGNs, the cosmological evolution of 
AGNs, and the distribution of the Extragalactic background light (EBL). 

The purpose of this paper is to consider the prospect of future blazar surveys by CTA. For this 
purpose, the blazar gamma-ray luminosity function (GLF) and spectral energy distribution (SED) 
are needed. The blazar GLF has been studied in detail in many papers (see e.g. [f|, |5]]). A new 
blazar GLFs [Q, |5|] has recently developed on the basis of the latest determination of the X-ray 
luminosity function of AGNs [Q, [7]] featuring a so-called luminosity density evolution. Another 
new feature of these GLFs is taking into account the blazar SED sequence. The blazar sequence 
is a feature seen in the mean SED of blazars that the synchrotron and inverse Compton (IC) peak 
photon energies decrease as the bolometric luminosity increases [@, g]. The key parameters in GLF 
have been carefully determined to match the observed flux and redshift distribution of EGRET 
blazars by a likelihood analysis. Predictions for the Fermi mission are in good agreement with the 
recent Fermi observed data [||, |To| ]. 

We utilize the model by []|] to predict the expected number and disuibutions of physical quan- 
tities of VHE blazars in a future CTA sky survey. We also discuss the expected high redshift blazar 
studies by CTA based on the first year Fermi source catalog. Throughout this paper, we adopt the 
standard cosmological parameters of {h,Q.M,0.\) = (0.7,0.3,0.7). 

2. VHE Blazar Survey with CTA 

There are various survey modes for VHE blazars such as a blank field sky survey and a follow- 
up survey for selected targets at other wavelengths. Here we mainly consider the case of a blank 
field sky survey mode following [|0. It enables us to find VHE sources without any biases from 
other wavelength observations. It should be noted that there is already a class of blazars which is not 
detected by Fermi but by current IACTs such as 1ES 0229+220 JTT|]. In the CTA era, the number 
of such sources is expected to increase. This kind of sources will be a new key to understanding 
the VHE gamma-ray emission mechanism in AGN jets, since the deabsorbed spectrum appear 
unusually hard which poses challenges to conventional acceleration and emission models (see e.g. 

'http://www.mpp.mpg.de/%7Erwagner/sources/ 
http://tevcat.uchicago.edu/ 
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Figure 1: Expected cumulative source counts as a function of the integral gamma-ray flux of VHE blazars. 
The five panels correspond to different photon energies, as indicated in the panels. The green curves are 
in the case of blank field sky survey; the blue curves are in the case of following up Fermi blazars with 
a sensitivity limit of 3 x 10~ 9 photons cm~ 2 s _1 for photon flux above 100 MeV [|l4||. The intergalactic 
absorption by EBL is taken into account for the solid curves, but not in the dashed curves. The 5-(7 detection 
limits of CTA for the array configuration E and 50 hr observing time are also shown. The blue solid curve in 
the panel of 10 TeV is shifted upward artificially by a factor of 1.2 for the purpose of presentation, because 
the blue solid and green solid curves totally overlap with each other. The horizontal thin solid line is the total 
expected number of Fermi blazars with the Fermi sensitivity given above. 



[12]). Furthermore, they will make a constraint on the intergalactic magnetic field strength by the 
limits to the expected cascade flux in the GeV energy range [0|. Thus, it is important to carry out 
an unbiased survey with CTA. The case of following up Fermi blazars is also presented in [^] in 
detail. 

To evaluate the expected outcomes of a survey with CTA, we need a blazar SED model and a 
GLF model. We utilize the recent GLF model by [||] which takes into account blazar SED sequence 
[§, ||]. This GLF model is consistent with the very recent Fermi observations ^[I^]. Fig. |] shows 
the cumulative source count distributions in the entire sky predicted by []|] above the following five 
energies of 30 GeV, 100 GeV, 300 GeV, 1 TeV, and 10 TeV as indicated in panels. In this figure, we 
adopt the EBL absorption model by [15], since gamma-ray photons from the universe suffer from 
the intergalactic absorption effect. The expected number of blazars will not change significantly 
when we adopt the EBL model [ 16], while it will decrease from our prediction when we adopt the 
model JT7|]. The differences among EBL models are addressed in ^ in detail. We also show the 
expected counts in the case of following up the Fermi blazars. 

Since observation time is limited, Table [T] shows the expected source counts with 100 hours of 
total observing time for various survey areas of 40, 400, and 4000 deg 2 using the EBL model [|15|] . 
Since the Field-of-View (FoV) of Large Size Telescope of CTA will be ~ 20 deg 2 [TO], we assume 
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Figure 2: Cumulative redshift distribution of blazars above 30 GeV for each CTA array configuration sys- 
tems (as labeled in the figure), assuming 50 hr observing time for each object and 20 degree zenith angle 
observation. The horizontal solid line corresponds to an expected count equal to one. 

FoV of CTA to be 20 deg 2 . Sensitivities for various observational time per FoV are calculated by 
using internal CTA tools. CTA will be able to detect ~ 50 and ~ 160 blazars with 1 year and 10 
years blank field sky survey assuming 1000 hrs of observing time for 1 year, respectively. The 
results show that a wide and shallow sky survey is favored for a fixed survey time. A wide field 
survey is expected to lead to discoveries of serendipitous objects. 



Energy range 


50 hr/FoV 
40 deg 2 


5 hr/FoV 
400 deg 2 


0.5 hr/FoV 
4000 deg 2 


>30 GeV 


0.36 


1.6 


5.4 


>100 GeV 


0.35 


1.5 


5.2 


>300 GeV 


0.18 


0.83 


2.8 


>1 TeV 


0.06 


0.25 


0.9 



Table 1: Expected blazar source counts for 100 hr blank sky survey by CTA. 

3. Distant Blazar Study by CTA 

CTA will also probe a deeper VHE universe than current IACTs do. Fig. |2| shows the cumula- 
tive redshift distribution of detectable blazars in the entire sky for various CTA array configurations 
(see [|l|] for details) predicted using the blazar GLF [g]. CTA will find ~ 20 blazars above z = 1 and 
the expected highest redshift will be at z ~ 1 .4 using the EBL model fllrjj]. These blazars will enable 
us to study the cosmological evolution of VHE AGN populations and to find EBL constraints at 
high redshift. 

Since Fermi is currently observing the entire GeV sky, we also evaluate the detectability of 
Fermi blazars in the first year catalog [J3j] using the internal spectrum simulation tool and the EBL 
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Figure 3: Expected spectrum of 1FGL J1344. 2-1723 blazar which is located at z = 2.49, assuming the array 
configuration B, 50 hr observing time, and 20 degree zenith angle. The red curve shows the intrinsic spec- 
trum without the EBL absorption. The gray curve shows the spectrum of the Crab nebula for comparison. 
The blue data points show the expected spectral data by CTA including the EBL absorption effects. The 
expected significance is 9.03-C7. 



model [[T|]. Fig. | shows the highest redshift Fermi blazar object at z = 2.49, 1FGL J1344-1723, 
detectable by CTA with 50 hr observing time including the intergalactic absorption effect. The 
expected significance above 30 GeV is 9.03-a with array B. However, we assume a simple power- 
law extrapolation from the Fermi spectrum. Further detailed SED fitting studies are required to 
validate this detectability. For other sources at high redshift, there are 17 Fermi blazar candidates 
for CTA above z = 1 with 50 hr observing time for each sources. 



4. Discussion and Summary 

In this paper, we discussed the expected number of blazars and highest redshift accessible with 
future CTA observation based on a recent blazar GLF model. For a blank field sky survey, a wide 
and shallow search will enable us to carry out an efficient survey. CTA will detect ~ 50 and ~ 160 
blazars with 1 year and 10 years blank field sky survey, respectively. With CTA, we should be 
able to find a blazar at z ~ 1.4 (~ 20 blazars above z = 1). We also estimated the detectable Fermi 
blazars. By simply extrapolating the GeV spectrum of the Fermi blazars in the first catalog, CTA 
is expected to detect a blazar at z = 2.49. This will enable us to study the cosmological evolution 
of VHE blazars. 

Although we discussed only blazars in this paper, other classes of AGNs are also expected as 
potential sources for CTA, such as core emission from radio galaxies [18], kpc jet emission JI^], 



low luminous AGNs By observing various AGN populations and high redshift sources, CTA 
will help to unify AGN populations and study the cosmological evolution of the VHE universe. 
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